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A principal 'supernova neutrino challenge' is the computational difficulty of six-dimensional neutrino radia- 
tion hydrodynamics. The variety of resulting approximations has provoked a long history of uncertainty in the 
core-collapse supernova explosion mechanism, but recent work highlighting low- mode convection and a newly- 
recognized instability in spherical accretion shocks may signal (yet another) resolution. As part of its goal of 
elucidating the explosion mechanism, the Terascale Supernova Initiative is committed to meeting the full com- 
plexity of the computational challenge. The understanding of supernova neutrino emission gained in detailed 
simulations provides a potential basis for learning about two major remaining unknowns in neutrino flavor mix- 
ing: the value of the mixing angle #13, and distinguishing between "normal" and "inverted" mass hierarchies. 



1. CORE-COLLAPSE SUPERNOVAE 

Shortly after the discovery of the neutron in the 
early 1930s, Baade and Zwicky declared, "With 
all reserve we advance the view that supernovae 
represent the transitions from ordinary stars to 
neutron stars, which in their final stages consist 
of extremely closely packed neutrons" £Q- This 
turned out to be true, at least for some 'core- 
collapse' supernovae (those of Type Ib/Ic/II); a 
black hole is another possible outcome. The 
dominant fleshing-out of the core collapse pro- 
cess in the last two decades has been the delayed 
neutrino-driven explosion mechanism |2I3| . 

A core-collapse supernova results from the evo- 
lution of a massive star. For most of their exis- 
tence, stars burn hydrogen into helium. In stars 
at least eight times as massive as the Sun (8 M ), 
temperatures and densities become sufficiently 
high to burn through carbon to oxygen, neon, and 
magnesium; in stars of at least ~ 10 Mq, burning 
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continues through silicon to iron group elements. 
The iron group nuclei are the most tightly bound, 
and here burning in the core ceases. 

The iron core — supported by electron degen- 
eracy pressure instead of gas thermal pressure, 
because of cooling by neutrino emission from car- 
bon burning onwards — eventually becomes unsta- 
ble. Its inner portion undergoes homologous col- 
lapse (velocity proportional to radius), and the 
outer portion collapses supersonically. Electron 
capture on nuclei is one instability leading to col- 
lapse, and this process continues throughout col- 
lapse, producing neutrinos. These neutrinos es- 
cape freely until densities in the collapsing core 
become so high that even neutrinos are trapped. 

Collapse is halted soon after the matter exceeds 
nuclear density; at this point ("bounce"), a shock 
wave forms at the boundary between the homol- 
ogous and supersonically collapsing regions. The 
shock begins to move out, but after the shock 
passes some distance beyond the surface of the 
newly-born neutron star, it stalls as energy is lost 
to neutrino emission and endothermic dissocia- 
tion of heavy nuclei falling through the shock. 

It is natural to consider neutrino heating as a 
mechanism for shock revival, because neutrinos 



1 



2 



dominate the energetics of the post-bounce evo- 
lution. Initially, the nascent neutron star is a 
hot thermal bath of dense nuclear matter, elec- 
tron/positron pairs, photons, and neutrinos, con- 
taining most of the gravitational potential energy 
released during core collapse. Neutrinos, hav- 
ing the weakest interactions, are the most effi- 
cient means of cooling; they diffuse outward on a 
time scale of seconds, and eventually escape with 
about 99% of the released gravitational energy. 

Because neutrinos dominate the energetics of 
the system, a detailed understanding of their evo- 
lution will be integral to definitive accounts of 
the supernova process. If we want to understand 
the explosion — which accounts for only about 1% 
of the energy budget of the system — we should 
carefully account for the neutrinos' much larger 
contribution to the energy budget. 

What sort of computation is needed to follow 
the neutrinos' evolution? Deep inside the newly- 
born neutron star, the neutrinos and the fluid 
are tightly coupled (nearly in equilibrium); but 
as neutrinos are transported from inside the neu- 
tron star, they go from a nearly isotropic dif- 
fusive regime to strongly forward-peaked free- 
streaming. Heating behind the shock occurs pre- 
cisely in this transition region, and modeling this 
process accurately requires tracking both the en- 
ergy and angle dependence of the neutrino distri- 
bution functions at every point in space. 

A full treatment of this six-dimensional neu- 
trino radiation hydrodynamics problem is a major 
challenge, too costly for contemporary computa- 
tional resources. While much has been learned 
over the years through simulation of model sys- 
tems of reduced dimensionality, there is as yet 
no robust confirmation of the delayed neutrino- 
driven scenario described above (see Sec. 0). 

Of interest are recent detections of a handful of 
unusually energetic Type Ib/c supernovae (often 
called 'hypernovae') in connection with gamma- 
ray bursts. Some show evidence for jets, and 
some do not. Determining the relative rates of 
jet-like hypernovae, non-jet-like hypernovae, and 
'normal' supernovae, and the possible associated 
variety of mechanisms, are important challenges. 

In summary, the details of how the stalled 
shock is revived sufficiently to continue plowing 



through the outer layers of the progenitor star 
are unclear. In normal supernovae, it may well 
be that some combination of neutrino heating be- 
hind the shock, convection, and instability of the 
spherical accretion shock leads to the explosion 
(see Sec. It is tempting to think that rotation 
and magnetic fields may play a more significant 
role in the rare jet-like hypernovae, perhaps giv- 
ing birth to 'magnetars,' the class of neutron stars 
with unusually large magnetic fields. 

2. HISTORY OF NEUTRINO RADIA- 
TION HYDRODYNAMICS 

From the description of the core-collapse super- 
nova process in the previous section, several key 
aspects of physics that a simulation must address 
can be identified: gravity, fluid dynamics (includ- 
ing equation of state and composition), neutrino 
transport and interactions, and perhaps magnetic 
fields. Ideally each of these aspects would be im- 
plemented in general relativity. 

Here the last two decades' progress on just one 
of these simulation aspects is sketched: the high 
dimensionality (three space and three momen- 
tum space dimensions — not to mention time de- 
pendence) of neutrino radiation hydrodynamics 
(see Table 1). The development of this aspect 
of the simulations is intertwined with important 
advances in the field, but of course does not repre- 
sent every insight relevant to the explosion mech- 
anism obtained via simulation or otherwise. 

I pick up the story in 1982, when simulations 
showing the stalled shock reenergized by neutrino 
heating on a time scale of hundreds of millisec- 
onds were first performed [2J- This was initially 
achieved in a simulation with a total of 2 dimen- 
sions (spherical symmetry, and energy-dependent 
neutrino transport). But with the introduction 
of full general relativity and a correction in an 
outer boundary condition jllj . it became clear 
that these models would not explode without a 
mock-up of a doubly-diffusive fluid instability in 
the newly-born neutron star that serves to boost 
neutrino luminosities |llll2ll8H4j — a simulation 
of effective total dimensionality "2.5" (see Ta- 
ble 1). That the necessary conditions exist for 
this particular instability to operate has been 
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Table 1 

Selected neutrino radiation hydrodynamics milestones in stellar collapse simulations studying the long- 
term fate of the shock. The "Yes" entries in the "Explosion" column are all marked with an asterisk as a 
reminder that questions about the simulations — described in the main text — have prevented a consensus 
about the explosion mechanism. "Total dimensions" is the average of "Fluid space dimensions" and u v 
space dimensions," added to V momentum dimensions." The abbreviation "N" stands for 'Newtonian,' 
while "PN" — for 'Post-Newtonian' — stands for some attempt at inclusion of general relativistic effects, 
and "GR" denotes full relativity. A space dimensionality in quotes — like "1.5" — denotes an attempt 
at modeling higher dimensional effects within the context of a lower dimensional simulation. For the 
fluid, this is a mixing- length prescription in the neutron star ("NS") or the heating region ("HR") 
behind the stalled shock. For neutrino transport, it indicates one of two approaches: multidimensional 
diffusion in regions with strong radiation/fluid coupling, matched with a spherically symmetric 'light 
bulb' approximation in weakly coupled regions ("thick/thin"); or the (mostly) independent application 
of a spherically symmetric formalism/algorithm to separate spatial angle bins ( "ray- by-ray" ) . 
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disputed |18I39I40| : and though related phenom- 
ena may operate |4Uj . more recent simulations 
with energy-dependent neutrino transport and 
true two-dimensional fluid dynamics indicate that 
fluid motions are either suppressed by neutrino 
transport 19 or have little effect on neutrino lu- 
minosities and supernova dynamics |3fil38| . 

Recognizing that the profiles obtained in spher- 
ical symmetry implied convective instabilities, 
and that observations of supernova 1987A also 
pointed to asphericities, several groups explored 
fluid motions in two spatial dimensions in the su- 
pernova environment in the 1990s. 

One class of simplifications allowed for neutrino 
transport in "1.5" or 2 spatial dimensions, but 
with neutrino energy and angle dependence in- 
tegrated out, reducing a five dimensional prob- 
lem to "1.75" or 2 effective total dimensions (see 
Table 1) |14I15I16| . These simulations exhibited 
explosions, and elucidated an undeniably impor- 
tant physical effect: a negative entropy gradi- 
ent behind the stalled shock results in convection 
that increases the efficiency of heating by neutri- 
nos. However, in the scheme of Table 1, the in- 
ability to track the neutrino energy dependence 
in these simulations could be viewed as a minor 
step backwards in effective total dimensionality. 
The energy dependence of neutrino interactions 
has the important effect of enhancing core delep- 
tonization, which makes explosions more difficult 
|7I41I42| ; this raised the question of whether the 
exploding models of the early- and mid-1990s 
were too optimistic. 

This concern about the lack of neutrino en- 
ergy dependence received some support from a 
simulation in the late 1990s involving a different 
simplification of neutrino transport: the impo- 
sition of energy-dependent neutrino distributions 
from spherically symmetric simulations onto fluid 
dynamics in two spatial dimensions |2()j . Unlike 
the simulations discussed above, these did not ex- 
plode, casting doubt upon claims that convection- 
aided neutrino heating constituted a robust ex- 
plosion mechanism. 

The nagging qualitative difference between spa- 
tially multidimensional simulations with different 
neutrino transport approximations motivated in- 
terest in the possible importance of even more 



complete neutrino transport: Might the retention 
of both the energy and angle dependence of the 
neutrino distributions improve the chances of ex- 
plosion, as preliminary "snapshot" studies sug- 
gested 43 30 ? Of necessity, the first such sim- 
ulations were performed in spherical symmetry, 
which nevertheless represented an advance to a 
total dimensionality of 3 (see Table 1). Results 
from three different groups are in accord: Spher- 
ically symmetric models of iron core collapse do 
not explode, even with solid neutrino transport 
|21I29I31| and general relativity |32l34j . Recently, 
however, it has been shown that the more mod- 
est oxygen/neon/magnesium cores of the lightest 
stars to undergo core collapse (8-10 M Q ) may ex- 
plode in spherical symmetry |23I24| . 

The current state of the art in neutrino trans- 
port in supernova simulations determining the 
long-term fate of the shock has been achieved by 
one of the above groups (centered at the Max 
Planck Institute for Astrophysics in Garching), 
who deployed their spherically symmetric energy- 
and angle-dependent neutrino transport capabil- 
ity |22j along separate radial rays, with partial 
coupling between rays |37| . Initial results — from 
axisymmetric simulations with a restricted angu- 
lar domain — were negative with regards to explo- 
sions (in spite of the salutary effects of convec- 
tion, and also rotation) [36] . apparently support- 
ing the results of Ref. [201 - A n explosion was 
seen in one simulation |44| in which certain terms 
in the neutrino transport equation corresponding 
to Doppler shifts and angular aberration due to 
fluid motion were dropped; this simulation also 
yielded a neutron star mass and nucleosynthctic 
consequences in better agreement with observa- 
tions than the "successful" explosion simulations 
of the 1990s |15I16| . arguably because of more ac- 
curate neutrino transport in the case of both ob- 
servables. The continuing lesson is that getting 
the details of the neutrino transport right makes 
a difference. 

In addition to accurate neutrino transport, low- 
mode (£ = 1,2) instabilities that can develop only 
in simulations allowing the full range of polar an- 
gles may make a subtle but decisive difference, as 
in an explosion recently reported by the Garch- 
ing group |38I24| . This achievement was presaged 
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by earlier studies demonstrating the tendency for 
convective cells to merge to the lowest order al- 
lowed by the spatial domain 05] and recognizing 
a new spherical accretion shock instability |46| 
(discovered independently in a different context 
in Ref . 07J ) . These global asymmetries may even 
be sufficient to account for observed aspherici- 
ties that have often been attributed to rotation 
and/or magnetic fields. 

Surely every 'Yes' entry in the explosion col- 
umn of Table 1 has been hailed in its time as 
'the answer' (at least by some!), and as a com- 
munity we cannot help hoping once again that 
these recent developments mark the turning of a 
corner; but important work remains to verify if 
this is the case. Several groups are committed 
to further efforts. For example, the Terascale Su- 
pernova Initiative (TSI, which includes authors of 
Ref. comprises efforts aimed at "ray-by-ray" 
simulations 0S| like those of the Garching group, 
as well as full spatially multdimensional neutrino 
transport, both with energy dependence only |49j 
and with energy and angle dependence • De- 
lineation of the possible roles of rotation and mag- 
netic fields are also being pursued by TSI. At least 
one other group is pursuing full spatially multidi- 
mensional neutrino transport |51I52| . 

3. SUPERNOVA NEUTRINOS AND 
FLAVOR MIXING 

As one who is engaged in the long-term devel- 
opment of a new code to be used in supernova 
simulations — with results from this code not yet 
in hand — I have spent this opportunity discussing 
our field's response to the challenge of neutrino 
radiation hydrodynamics in core-collapse super- 
nova simulations over the last two decades, to- 
gether with many references. Hopefully this pro- 
vides some understanding of the history of claims 
in our field, and gives an appreciation of the years 
of dedicated work by many individuals to advance 
the state of the art. But in closing, I should at 
least say something more directly tied to the topic 
of the conference! 

So how do these large-scale simulations relate 
to neutrino flavor mixing? At present, the sim- 
ulations assume massless neutrinos, relying on 



the standard expectation that the effective mass 
from neutrino forward scattering off electrons 
suppresses flavor mixing in the high-density re- 
gion where the explosion is launched. 2 Before 
reaching terrestrial detectors, however, neutri- 
nos will pass through flavor-changing resonances 
at lower densities in the stellar envelope (and 
possibly also experience flavor-changing effects 
while traversing Earth before entering a detec- 
tor) . Hence collapse simulations provide neutrino 
spectra as input to studies that take these effects 
into account in determining possible signals in 
terrestrial neutrino detectors (see, for example, 
overviews in Refs. |53l54j ). To the extent that 
generic energy- and time-dependent features of 
the signals can be motivated by the large-scale 
simulations, such studies provide a possible basis 
for learning — from neutrino detection of a future 
Galactic supernova — about two major remaining 
unknowns surrounding neutrino flavor mixing: 
the value of the mixing angle #13, and distinguish- 
ing between "normal" and "inverted" mass hier- 
archies. Alternatively, to the extent that these 
neutrino mixing features are determined by ter- 
restrial long-baseline neutrino experiments before 
the next Galactic supernova, neutrino signals will 
allow the details of supernova dynamics to be ob- 
served in real time. 
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